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Abstract 

In our recent study we have determined two new spatially dependent nuclear PDF (nPDF) sets, 
EPS09s and EKS98s. With these, the hard-process cross-sections can be calculated in different 
centrality classes consistently with the globally analyzed nPDFs for the first time. The sets were 
determined by exploiting the A-systematics of the globally fitted nPDF sets, EPS09 and EKS98. 
For the spatial dependence of the nPDFs we used a power series ansatz in the nuclear thickness 
function Ta- In this flash talk we introduce the framework, and present our NLO EPS09s-based 
predictions for the nuclear modification factor in four centrality classes for inclusive neutral pion 
production in p+Pb collisions at the LHC and for inclusive prompt photon production in d+Au 
collisions at RHIC at midrapidity. 



1. Introduction 

According to the collinear factorization theorem |Q] [2J the inclusive yield of a parton A; in a 
heavy ion collision at a given impact parameter b can be calculated from 

dN AB^k + x = TM ff( Xl ,Q 2 ) ® ff( X2 , Q 2 ) ® d«x<'^ +r + 0(l/Q 2 ), (1) 

i,j,X' 

where T AB (b) is the nuclear overlap function, ff{x\,Q 2 ) (/^(*2>2 2 )) is a universal, process- 
independent nuclear parton distribution function (nPDF) describing the momentum distribution 
of a parton i (j) inside the nucleus A (B), and d&'^ k+x can be computed using perturbative 
QCD. The nPDFs are known to be modified with respect to corresponding free proton PDFs and 
this modification can be quantified by defining the ratio Rf(x, Q 2 ) as 

fif(i,e 2 )E^(i,e 2 )/^(i,e 2 ). (2) 

As the (n)PDFs are of non-perturbative origin, the nuclear modifications of the PDFs cannot be 
computed from the first principles of QCD, but these can be determined by global fitting to the 
experimental data using DGLAP evolution as done e.g. in Refs. J[3j H) . So far all the global fits 
have considered only minimum bias observables and thus the resulting nPDFs should be con- 
sidered as spatially averaged quantities. This implies that one has not been able to consistently 
compute the hard-process cross sections in different centrality classes. In [5| we have considered 
this problem and developed a framework for spatial dependence of the nPDFs. 
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2. Framework 



We introduce a nuclear modification rf(x, Q 2 , s) which now depends also on the transverse 
position s in the nucleus. For this we require that the spatial average is identical to the globally 
analyzed nuclear modification 

«?(*, Q 2 ) = \ J d 2 s T A (s) rf(x, Q\ s), (3) 

which we take here from the earlier global analysis, EPS09 [3 1 or EKS98 [4|. As this constrains 
only the spatially averaged modification, we need to make an assumption about the spatial depen- 
dence. Motivated by the small-x shadowing studies @, we assume that the nuclear modification 
is related to the nuclear thickness function Ta(s) and choose a power series form 

n 

rf{x, Q\ s) = 1 + £ c)(x, Q 2 ) [T A (s)V, (4) 
;=i 

where c'j{x, Q 2 )'s are now A-independent fit parameters. We determine the values of these pa- 
rameters by requiring the same A-dependence for the spatially averaged modification as obtained 
in the earlier global analysis, i.e. we minimize the^ 2 defined as 

n 2 



*?(*. Q 2 ) = 2 



Rf(x, Q 2 )-{J d 2 s 7U(s) rf{x, Q 2 , s) 



Wf(x, Q 2 ) 



(5) 



For the weight factor we use Wf (x, Q 2 ) = 1 for EPS09 and Wf(x, Q 2 ) = 1 - Rf(x, Q 2 ) for 
EKS98. In figure [T] we plot the spatially averaged gluon modification as a function of A from 
EPS09LO and EPS09NLO and from the corresponding spatial fits for fixed x and Q 2 values. We 
observe that taking n — 4 in equation Q is sufficient to reproduce the A-systematics accurately 
for all the sets considered. By repeating this fitting procedure over all parton flavors and over 
an (x, Q 2 ) grid for all different sets (EKS98: LO, EPS09: LO and NLO + 30 error sets each) 
we obtain the new spatially dependent nPDF sets EPS09s and EKS98s, which are now available 
at our webpag^] In figure [2] we plot the gluon modification from EPS09sNLO as a function 
of x and s{— |s|) at the EPS09 initial scale Q 2 = 1.69 GeV 2 . From these two figures one can 
notice how the A-dependence of the Rf(x, Q 2 ) maps to the spatial depedence of the rf(x, Q 2 , s): 
modifications are larger at the thick center of nucleus and decrease towards the thin edge. 



3. Applications 

In Ref. we compared our NLO calculation of the nuclear modification factor for n° pro- 
duction in d+Au collisions at y — and V^wv = 200 GeV with the PHENIX measurement 
|7| in four centrality classes. Within the uncertainties, the calculations described the data well 
in all these centrality classes. In figure [3] we plot similar calculation for p+Pb collisions at 
s/snm = 5.0 TeV and y — in nucleon-nucleon CMS-frame at the LHC. The centrality classes 
we consider here are - 20%, 20 - 40%, 40 - 60%, and 60 - 80%. These are determined 
in terms of impact parameter intervals which we calculate using the optical Glauber model, see 
Ref. for details. The calculations are done with the INCNLO-packageJ^] HI and CTEQ6M |9| 
free proton PDFs with the EPS09s modifications, and with three fragmentation functions (FFs): 
KKP lITOl . AKK ifTTI . and fDSS lfT2l . Our result for the corresponding minimum-bias ratio 
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Figure 1: The spatially averaged gluon modification 
R g (x, Q 2 ) from EPS09 and EPS09s in LO and NLO as a 
function of A for fixed values of x and Q 2 . From 1 5 1 . 



Figure 2: The spatially dependent gluon modification 
r g (x, Q 2 ,s) from EPS09sNLO for A = 208 as a function 
of x and s at the initial scale of EPS09. From |5 1. 



R n ° vh (p T ) (left panel of figure 17 in Ref. Q) agrees encouragingly well with the first ALICE data 



for Rf pb (p T ) ED. From figure^ we observe that the nuclear effects are more pronounced in the 
most central collisions and they decrease towards more peripheral collisions. Also, we see that 
the possible differences between different FFs essentially cancel out in ratios like R^ pb (pT). 

In figure |4] we plot the nuclear modification factor for inclusive prompt (= direct + fragmen- 
tation) photon production in d+Au collisions at -\Js NN = 200 GeV and y = 0. For the parton 
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Figure 3: The nuclear modification factor i?^ pb for inclusive neutral pion production in p+Pb collisions at the LHC for 
four different centrality classes (different panels) at midrapidity, calculated with three different FFs, KKP (dot-dashed), 
AKK (dashed) and fDSS (solid). The blue uncertainty band is calculated from the 30 error sets in EPS09s with fDSS 
FFs. From |5 1. 
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fragmentation to photons we have used the BFG (set II) FFs [14|. To quantify the isospin ef- 
fect we have plotted also the ratio R^ Aa without the nuclear modifications to each panel. At 
Pt ~ 5 GeV/c the antishadowing in the nPDFs seems to compensate for the isospin suppression 
and at pj ~ 15 — 20GeV/c the isospin suppression dominates so the centrality dependence of 
R y ,. is minimal at this region. 
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Figure 4: The nuclear modification factor R y iAll for inclusive prompt photon production in d+Au collisions at RHIC for 
four centrality classes (different panels) at midrapidity with and without the nuclear modification of the PDFs. Frag- 
mentation component is calculated with BFG FFs and the blue uncertainty band is calculated from the 30 error sets in 
EPS09s. From fEl . 
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